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We study the reduced single-particle density matrix (RSPDM), the momentum distribution, nat- 
ural orbitals and their occupancies, of dark "soliton" (DS) states in a Tonks-Girardeau gas. DS 
states are specially tailored excited many-body eigenstates, which have a dark solitonic notch in 
their single-particle density. The momentum distribution of DS states has a characteristic shape 
with two sharp spikes. We find that the two spikes arise due to the high degree of correlation 
observed within the RSPDM between the mirror points (x and —x) with respect to the dark notch 
at x — 0; the correlations oscillate rather than decay as the points x and —x are being separated. 

PACS numbers: 03.75.-b,03.75.Kk 
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I. INTRODUCTION 

Exactly solvable models have the possibility of pro- 
viding important insight into the quantum many-body 
physics beyond various approximation schemes. Two 
of such models, the Tonks-Girardeau [l[ and the Lieb- 
Liniger model Q, which describe interacting Bose gases 
in one dimension (ID), have drawn considerable attention 
in recent years with the developments of the experimen- 
tal techniques for tightly confining atoms in effectively ID 
atomic waveguides [H, 0, [f| |(| . The Lieb-Liniger model 
(LL) describes a system of bosons interacting via two- 
body (5-function interactions 0. The Tonks-Girardeau 
(TG) model corresponds to infinitely repulsive ("impen- 
etrable core") bosons in ID [3, 0); this model is exactly 
solvable via Fermi-Bose mapping, which relates the TG 
gas to a system of noninteracting spinless fermions in 
ID A study of atomic scattering for atoms confined 
transversally in an atomic waveguide has lead to a sug- 
gestion for the experimental observation of a TG gas [8| ; 
such atomic systems enter the TG regime at low tem- 
peratures, low linear densities, and strong effective inter- 
actions [a [t| [l(|. The experimental realization of the 
TG model was reported in two experiments from 2004 
0,0]. Moreover, nonequilibrium dynamics of the ID in- 
teracting Bose gases including the TG regime has been 
recently experimentally studied within the context of re- 
laxation to an equilibrium Q . Within this paper we ana- 
lyze the reduced single-particle density matrix (RSPDM) 
and related observables of certain specially tailored ex- 
cited eigenstates of the TG gas, which are also referred 
to as dark "soliton" (DS) states [II [II [II. 

Dark solitons are fundamental nonlinear excitations. 
Within the context of interacting Bose gases, they were 
mainly studied in the regime of weak repulsive inter- 
actions [3, [H, [H, [TtJ were mean-field theories [e.g., 
the Gross-Pitaevskii theory, which employs the nonlin- 
ear Schrodinger equation (NLSE)] are applicable. In the 
regime of strong repulsive interactions in quasi-lD geom- 
etry, dark solitons were also studied by using NLSE with 
a quintic nonlinear term [H, QjJ [2(| • In Ref. [H| , Gi- 



rardeau and Wright have studied the concept of dark soli- 
tons within the exactly solvable TG model; they found 
specially tailored excited many-body eigenstates of the 
TG gas on the ring (DS states), with a dark notch in 
their single-particle density, which is similar to the dark 
notch of nonlinear dark-solitons. The dynamics of such 
excitations in a TG gas was studied by Busch and Huyet 
[12| in a harmonic trap. Recently, a scheme based on 
parity selective filtering (" eva por ation" ) of a many-body 
wave function was suggested [13J as a candidate for the 
experimental observation of DS states. However, to the 
best of our knowledge, the momentum distribution, the 
RSPDM, natural orbitals (NOs) and their occupancies, 
have not been studied yet for DS states. These quan- 
tities are important for the better understanding of DS 
states, but may also be necessary ingredients for their ex- 
perimental detection, which provides motivation for this 
study. 

The calculation of correlation functions (such as the 
RSPDM) 




27, 28, 29, 301, M, |3J, |33j, |3J, |3£ 
40, 41 , [42 , [43| from the many-body wave functions 
|U Is Hi li 144 45, Hil yields important physical informa- 
tion (such as the momentum distribution) on the state of 
the system. Within the TG model, the RSPDM and the 
momentum distribution have been studied in the con- 
tinuous [SI, [H HIj27|, JM M M M El, and dis- 
crete (lattice) case [31], |34|, |33 S3, both for the 
static 




dependent problems 



381 . [43|| and time- 
In the stationary 



case, most studies consider the ground state properties of 
the TG gas. The momentum distribution for the ground 
state of the TG gas on the ring has a spike at k = 0, 
nn(k) oc |fc| -1 / 2 [2l|. In both the harmonic confinement 
[ij, [2^| and on the ring [26J , the TG ground state mo- 
mentum distribution decays as a power law ns(Ai) cx fc~ 4 ; 
in Ref. [26| it has been pointed out that fc _4 -decay is also 
valid for the LL gas (for any strength of the interaction) . 
These grou nd states of the TG gas are not Bose con- 
densed [HI, [28[, which is evident from the fact that the 
occupancy of the leading natural orbital scales as \/~N for 
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large N [27, 28]. In the box-confinement, the momentum 
distribution of a TG gas has been studied by generalizing 
the Haldane's harmonic-fluid approach [25|]. Besides for 
the ground states, the momentum distribution has been 
analyzed in time-dependent problems including irregu- 
lar dynamics on the ring [33l |. dynamics in the harmonic 
potential with time dependent frequency 35], and in a 
periodic potential in the context of many-body Bragg 
reflections [40] . A number of interesting results for time- 
dependent problems have been recently obtained within 
the discrete lattice model including fermionization of the 
momentum distribution during ID free expansion [34| . 
and relaxation to a steady state carrying memory of ini- 
tial conditions [39]. 

The correlation functions for TG and LL models were 
studied by using various a naly tical and numerical meth- 



ods 21 



22, 23 



24J25 



2G 



H, 11, [33, H, M, SO, \M El IIILThe formula that was 
derived and employed in Ref. [40l | allows efficient and 
exact numerical calculation of the RSPDM for the TG 
gas in versatile state s (g round state, excited eigenstates, 
time-evolving states [40] ) , and for a fairly large number of 
particles. We find it suitable for this study of DS states. 

Here we numerically calculate the RSPDM correla- 
tions, natural orbitals and their occupancies, and the mo- 
mentum distribution of DS states. We find that these ex- 
cited eigenstates of a TG gas have characteristic shape of 
the momentum distribution with two sharp spikes. The 
two sharp spikes arise due to the high degree of corre- 
lation observed within the RSPDM between the mirror 
points, x and —x, with respect to the dark notch at x = 0; 
interestingly, the correlations oscillate rather than decay 
as the points x and —x are being separated. 



II. THE MODEL 

We study a system of N identical Bose particles in 
ID space, which experience an external potential V(x). 
The bosons interact with impenetrable pointlike inter- 
actions [l|, which is most conveniently represented as a 
subsidiary condition on the many-body wave function [l[ : 



ip B (xi,x 2 , ■ ■ ■ ,x N ,t) = if Xi = Xj 



(1) 



for any i ^= j. Besides this condition, ips obeys the 
Schrodinger equation 



' dt 



N 

£ 



3? 

dx 



2 +V( Xj ) 



(2) 



here we use dimensionless units as in Ref. [13j, i.e., 
x = X/X , t = T/T , and V(x) = U(X)/E , where 
X and T are space and time variables in physical units, 
Xq is an arbitrary spatial length-scale (e.g., Xq — 1 /im), 
which sets the time-scale Tq = 2mX$ j h, and energy-scale 



Eq = h 2 / (2itiXq); m denotes particle mass, and U(X) is 
the potential in physical units. The wave functions are 
normalized as J dx\ . . . dxN\4>B{xi,x 2 , ■ ■ ■ , xjy, t)\ 2 = 1< 

The solution of this system may be written in compact 
form via the famous Fermi-Bose mapping, which relates 
the TG bosonic wave function ips to an antisymmetric 
many-body wave function ipp describing a system of non- 
interacting spinless fermions in ID [l|: 



V>b = A(xi, . . . ,x N )ip F (xi,x 2 , ■ ■ ■ ,x N ,t). (3) 



Here 



A = ni<i<j<Arsgn(xi 



(4) 



is a " unit antisymmetric function" [1| , which ensures that 
ipB h a s proper bosonic symmetry under the exchange of 
two bosons. The fermionic wave function ipp is com- 
pactly written in a form of the Slater determinant, 



iPf(xi, ■ ■ ■ ,%N,t) 



1 



JY 



det [ip m (xj,t)], (5) 



where ip m (x,t) denote N orthonormal SP wave func- 
tions obeying a set of uncoupled single-particle (SP) 
Schrodinger equations 



.dipm 



d 2 



ip m (x,t), m = 1, ... ,7V. (6) 



Equations ©-© prescribe construction of the many- 
body wave function describing the TG gas in an external 
potential V(x), both in the static [l| and time-dependent 
case [TTj l . The eigenstates of the TG system are 



i/>b(xi, . . - ,x N ) = A(xx 



1 N 

,XN)-r= det IfimiXj)}, 

V-/V! m J =1 



(7) 

where (f> m (x) are single-particle eigenstates for the po- 
tential V(x). In the rest of the paper we will discuss 
the eigenstates of the TG system and their observables; 
hence, we drop the time-variable from subsequent nota- 
tion. 

The many-body wave function %Pb fully describes the 
state of the system. However, its form does not transpar- 
ently yield physical information related to many impor- 
tant observables (e.g., the momentum distribution). The 
expectation values of one-body observables are readily 
obtained from the RSPDM, defined as 



PB{x,y) = N jdx 2 ■ ■ .dx N i> B {x,x 2 , ■ ■ ■ ,x N )* 

xip B (y,x 2 ,...,x N ). (8) 
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The observables of great interest are the SP x-density 
Pr(x, x) = m —i \<j>m(%)\ 2 i an d the momentum distri- 
bution [211 ]: 

n B (k) = ^- Jdxdye lk( - x -y'>p B (x,y). (9) 

The SP density pb(x,x) is identical for the TG gas and 
the noninteracting Fermi gas [l[, however, the momen- 
tum distributions of the two systems considerably differ 

a 

A concept that is very useful for the understanding of 
the many-body systems is that of natural orbitals (NOs) . 
The NOs are eigenfunctions of the RSPDM, 

[dxp B (x,y)<5>i(x) = A 4 $,(y), i = 1,2,..., (10) 

where Xi are the corresponding eigenvalues; the RSPDM 
is diagonal in the basis of NOs, 

oo 

PBfotf) =£*i*i (*)*<(»)• (11) 

The NOs can be interpreted as effective SP states occu- 
pied by the bosons, where represents the occupancy 
of the corresponding NO [HI]. The sum of the Fourier 
power spectra of the NOs is the momentum distribution: 



n fl (fc) = ^Ajij(fc)i 4 (fc), (12) 

i=l 

where is the Fourier transform of 3>i(x). 

The RSPDM of the noninteracting fermionic system 
on the Fermi side of the mapping is 

N 

pf(x,v) = ^*m{ x )4>m{y); (13) 

m— 1 

evidently, the SP eigenstates <fi m (xj) are NOs of the 
fermionic system, with occupancy unity [23|. The 
fermionic momentum distribution is 

N 

n F (k) = J2 &,(*0M*). (14) 

m— 1 

where 4> m (k) is the Fourier transform of <f> m (x). 

The calculation of the TG momentum distribution is 
preceded by a calculation of pb{x,u), which we conduct 
according to the method described in Ref. [4(|. If the 
RSPDM is expressed in terms of the SP eigenstates <p m 
as 



N 

Pb{x,v) = 4>*i{x)A ij [x,y)(l) j {y), (15) 

it can be shown that the N x N matrix A(x, y) — 
{Aij(x,y)} has the form 

A(x,y) = (P~ 1 ) T detP, (16) 

where the entries of the matrix P are Py (x, y) = <5y — 
2 (J dx'4>*(x')4>i(x') (x < y without loss of generality) 
pfjf . Formulas (fTS"]) and (|16p enable fast numerical cal- 
culation of the RSPDM (and related quantities) for dark 
"soliton" states. 



III. DS STATES ON THE RING 

Within this section we analyze the RSPDM, the mo- 
mentum distribution, NOs and their occupancies for ex- 
cited eigenstates of a TG gas on the ring of length L; 
in other words, external potential is zero, x-space is 
x € [— L/2, L/2], and periodic boundary conditions are 
imposed. The many-body eigenstates of the TG gas 
are constructed from the SP eigenstates of the system 
via Eq. (J7J). Th SP eigenstates for the ring geometry 
are plane waves \ / Le lkmX , with SP energy k^; here 
k m = 2jTm/L, and m is integer (48j . Apparently, the 
eigenstates yjl/ Le lkmX and 1 / Le~ lkmX are degener- 
ate. This degeneracy in the SP eigenstates induces [via 
Eq. 0] degeneracy of the TG many-body excited eigen- 
states. One particular subspace of degenerate eigenstates 
(DEs) is spanned with 

4> m {x) = -L[a-e-* k ™ x + a+e^ x ] (17) 

where |a~ | 2 + |a+ | 2 = 1, and m = 1, . . . , N; the corre- 
sponding many-body eigenstates are 

N 

1pDB = A(xi, . . . ,X N )L~~ X 

det [a-e- lkmX * + a+e lkmX ^}. (18) 

j,m=l 

Intuition suggests that, although these states are degen- 
erate, some of the corresponding observables, such as the 
SP density in x-space, the momentum distribution, spa- 
tial coherence or entropy, could be quite different from 
one eigenstate to another depending on their internal 
symmetry, which is designated by the choice of the coef- 
ficients a~ and a+. 

In Ref. [ll[, Girardeau and Wright have pointed out 
that if one constructs excited many-body eigenstates of 
the TG gas on the ring as 
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FIG. 1: (color online) The RSPDM of dark "soliton" states, 
PDs(x,y), for N = 5 (a), 11 (b), 17 (c), and 25 (d). 




FIG. 2: (color online) The momentum distributions corre- 
sponding to DS states (a), and to degenerate eigenstates with 
randomly chosen phases (b); figures are shown for N = 11 (x- 
symbols, blue dotted line), 17 (diamonds, red dashed line), 
and 25 (circles, solid black line). 



analytical form up to N = 7 by using Mathematica. How- 
ever, for larger N we resorted to numerical calculations. 
It is straightforward to see that the RSPDMs of two DS 
states, for two different values of L, say L\ and L 2 , are 
connected by a simple scaling, 



tpDS = A(xi, . . -,x N ) 



N 



det [sin k m Xj], 

3)711=1 



(19) 



that is, if one chooses the coefficients as a m — i/y/2 and 
a m = ~zi/v^2, the SP density of these many-body eigen- 
states irHl. 



PDS(X,X) 



N + 1 sin( 



(jV+ i )2 H cos(^ 



Lsin(2|2i) 



(20) 



has the structure closely resembling dark solitons 
(hence the notation iJjds for the many-body wave func- 
tion, and analogously for related observables below). The 
structure of these excited eigenstates is somewhat ar- 
tificial because on the fermionic side of the mapping, 
these states correspond to noninteracting fermions be- 
ing placed solely within the odd SP eigenstates sin kffiX . 
Nevertheless, such states can be excited by filtering of 
the many-body wave function [l3j . 

Let us utilize the procedure outlined in Sec. II to cal- 
culate the RSPDM, and related one-body observables for 
DS states [Eq. ([TO)) ]. It is straightforward to calculate 
the entries of the matrix P = 1 — Q [ see Eq. (|16p ] , where 



Qi 



sin(2(i + j)nx/L) sin(2(? — j)irx/L) 



sin(2(i + ])ny/L) 



(i + j)ir 



(i - j)ir 
sin(2(i - j)ity/L) 
(i - j)n 



Qii ' ■ 2 - 



(1H£) sin(^) 



L 



2iir 



2in 



(21) 



for i, j 



1. 



, N. As for the inverse of the matrix P, 



and consequently the RSPDM, we were able to find its 



LiPDS.LtixLi^yLi) = L 2 pDS,L 2 (xL2,yL 2 ), 



(22) 



where x,y € [— h , thus, it is sufficient to calculate it 
for just one value of L. In what follows, without loosing 
any generality, we choose N = L. 

Figure [1] displays contour plots of Pds( x jU) f° r N = 
5, 11, 17 and 25. We clearly see a characteristic pattern 
for each value of N: The RSPDMs are largest close to the 
diagonal, with oscillations following the x-space density 
from Eq. (f20|) . However, there are strong correlations 
along the line x — — y indicating coherence between mir- 
ror points x and —x around the DS center (at x = 0). 

Figure [2Ja) displays the momentum distribution 
n Ds(k) of DS states for N = 11, 17, and 25. All momen- 
tum distributions for the ring geometry are normalized 
as Sfe n B{k m ) — N (the SP momentum values k m are 
discrete in the ring geometry). The momentum distri- 
butions have a characteristic shape with a smooth hump 
close to the origin (k — 0), and with two sharp spikes 

which are located at ±k pea k — ± J2m=i km/N — ±7r(A+ 
1)/L; the spikes indicate that there is high probability of 
finding a boson in momentum states exp(±z7r(A r + l)/L). 
Note that due to our choice N = L the peaks for differ- 
ent values of N approximately coincide at ±7r(l + l/A) rj 

±7T. 

The sharp spikes at /c pea fc = ±7r(A + V)/L are in- 
timately related to the strong correlation between the 
mirror points x and —x. This is illustrated in Fig. [3] 
which shows the cross-diagonal section of the RSPDM 
Pds(x, —x) and the function cos(2/c peQJ tj:) for N = 25. 
There is evident correlation between pos{x,—x) and 
cos(2fc pea fea:). Because of the symmetry pus{x,y) — 
PDs(y,x), the Fourier transform (FT) with respect to 
exp[ik(x — y)] reduces to FT with respect to cos k(x — y), 
which is cos2fcir at y = —x; hence, from Fig. [3] it 
immediately follows that the cross-diagonal behavior of 
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FIG. 3: (color online) The cross-diagonal Pds{x,— x) of the 
dark "soliton" RSPDM (solid black line), displaying long- 
range oscillatory correlations between mirror points x and 
— x, and the cosine function cos(2k pea kx) (red dotted line). 



FIG. 5: (color online) The occupancies of the NOs for the 
state tfiDS (black diamonds) and a typical state ipDE (red 
squares), for N = 25 particles. The sharp drop in the occu- 
pancies of tpDS occurs between Ajv and Ajv+i- 
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FIG. 4: (color online) The RSPDM of typical eigenstates 
ipDE [Eq. (]23)] for N = 7 (a), 11 (b), 17 (c), and 25 (d). 



Pds{x,—x) induces the peaks in the momentum distri- 
bution of DS states. We would like to point out that 
the correlations Pds(x, ~ x) do not decay, but oscillate, 
as the separation between points x and —x is increased. 

In order to gain more insight into the origin of the 
two sharp spikes in the momentum distribution and the 
related coherence between mirror points x and —x, it 
is illustrative to calculate the RSPDM and the momen- 
tum distribution for eigenstates that are degenerate (i.e., 
that have the same energy) to DS states, but which are 
less restrictive with respect to symmetry of the coef- 
ficients a~ and a+. If the coefficients are chosen as 
a~ = i exp(i9 m ) / \/2 and a+ = — i exp(— id m )/y/2, one 
obtains a whole class of eigenstates degenerate to dark 
solitons, which have the form 



— 

ip DE = A(x x , . . . ,X N ) (^-J X 
N 

det [sin(k m Xj + 9 m )], (23) 

j,m= 1 

where 9 m , m — 1,...,N are N phases (for 9 m = 0, 
ipDS = i>DE)- 

Figure 2] displays contour plots of RSPDMs pde(x,v), 
which corresponds to some typical states ipDE obtained 
from Eq. ([23"]) by randomly choosing N phases 9 m (with 
respect to the uniform probability density in [— 7r,7r]). 
We see that the SP density for this state, pde(x,x), is 
not zero at x = 0, which evidently follows from the fact 
that shx{k m x + 9 m ) is not an odd function for 9 m ^ 0, 
while pde(x,x) = J2m=i I sin(fc m a;-|-6 l m )| 2 . Furthermore, 
we observe that the structure of the RSPDM along the 
x = —y line is absent, that is, there is no coherence 
between the mirror points x and —x. A closely related 
observation is that the momentum distributions of such 
states do not have a pair of sharp spikes which are present 
in riDs(k); this is illustrated in Fig. [2Jb) which shows 
typical momentum distributions fiEE(k) for N — 11, 17 
and 25. 

Besides the RSPDM and the momentum distribution, 
excited many-body eigenstates of interest can be charac- 
terized by the corresponding natural orbitals (NOs) and 
their occupancies. Figure [5] shows the occupancies of 
the NOs of the state ipDS, and a typical state ipDE, for 
N = 25. We observe that the occupancies are fairly low 
(less than one) for all NO's, but there is a sharp drop 
in the occupancies after the 25th NO. We have observed 
such a behavior for other values of N as well. In contrast, 
the occupancies of the NOs corresponding to a typical 
state ipDE do not exhibit a sharp drop after the A^th 
orbital, but decrease rather smoothly. 



6 



i (a) | 

, J 






\ 


-10 -5 X 5 10 




-10 k 10 


■ A 

(d) [ / \ 

A J 


i \ 



-10 -5 x 5 10 



FIG. 6: (color online) The NOs of a DS state in a;-space (left 
column) and their power spectra in fc-space (right column) 
for N = 25. (a) The first $\{x) (black solid line), and the 
third $z(x) (red dotted line) NO. (b) The Fourier power spec- 
trum of the first |$i(fc)| 2 (black circles, solid line), and the 
third |$3(fc)| 2 (red diamonds, dotted line) NO. (c) The 24th 
(red dotted line) and 25th (black solid line) NO in z-space. 
(d) The momentum distribution (red squares, dotted line) in 
comparison to the contribution from the 24th and 25th NO: 
5!^ i=24 Ai|<£ , i(fc)| 2 (black circles, solid line). 



Figure [5] illustrates the spatial structure and the 
Fourier power spectra of the NOs corresponding to the 
DS state for N = 25. The spatial structure of calcu- 
lated NOs is either symmetric or antisymmetric. This 
is connected to the symmetry pds(x, y) = Pds{~V, —x); 
due to this symmetry it follows that if some NO is non- 
degenerate, it is either symmetric or antisymmetric; if 
two NOs are degenerate (their occupancies are identi- 
cal), they can be superimposed to yield one symmetric 
and one antisymmetric NO. Our numerical study shows 
that the low order (leading) NOs are localized in space, 
but broad in fc-space; Fig. [5Ja) depicts the £-space struc- 
ture, and Fig. b) show the fc-space structure of the first 
and the third NO. We see that these low order NOs do 
not contribute to the sharp peaks observed in the mo- 
mentum distribution of DS states. Further inspection of 
the NOs reveals that those NOs just on the upper side of 
the sharp drop in \j (Fig. [5]) are in fact responsible for 
the sharp peaks: Figs. [S^c) and (d) display the x-space 
and fc-space structure, respectively, of the 24th and the 
25th NO (N = 25). The total momentum distribution 
(red squares, dotted line in Fig. [H£d)) can be written 
as YliZi a contribution to this sum stem- 

ming from the 24th and the 25th NO is shown in Fig. 
|6fd) with black solid line. Evidently, for this DS state 
where N = 25, the 24th and the 25th NO give rise to the 
peaks in the momentum distribution. 

It is interesting to note that when all phases are chosen 
to be identical but not zero, e.g., if 9 m — tt/2, then all 
of the fermionic NOs are cx cos(fc m x), we again observe a 
higher degree of correlation between mirror points in the 



N 


S[^D S ] 


S[lpDE] 


11 


2.90 


3.21 


17 


3.39 


3.66 


25 


3.83 


4.05 



TABLE I: The entropy S of dark "soliton" states, and typical 
tpnE states for different values of the number of particles N. 



RSPDM and peaks in the momentum distribution (not 
shown). 

All of the observations above indicate a somewhat 
smaller degree of order in the degenerate eigenstates i/jde 
than in dark solitons ipDS, which follows from the ran- 
dom (disordered) choice of the phases 9 m . This is fur- 
ther underpinned in Table [H which shows the entropy 
S = — J2iPi log where pi = Xi/N, for the dark "soli- 
ton" states ipDS, and typical ipDE states. The entropy 
of states ipDE is systematically larger than in the states 

"0-DS- 

From our observations it follows that the many-body 
state ipDS contains a distinct component, which can be 
interpreted as a standing wave populating momentum 
modes at ±k pea k- In the effective single-particle picture, 
we see that this component give rise to the population of 
the natural orbitals close to (and including) the Nth NO. 
However, it should be pointed out that this component 
is fairly small, i.e., it yields small occupation of these 
effective SP states. 

In a similar fashion to the excited ipDS state, the 
ground state of the TG gas on the ring yields distinct 
population of the zero- momentum mode [2l| ; in this case, 
however, the zero-momentum mode is the leading natu- 
ral orbital, and its population is fairly large (it scales as 
VN [HI). Even though the TG states are not Bose con- 
densed, they can sharply populate a single momentum 
mode because bosons do not obey the Pauli principle and 
consequently more than one boson can occupy a single 
momentum state (which is not the case for noninteracting 
fermions). 

It is interesting to note that on the Fermi side of the 
mapping, the momentum distribution of noninteracting 
fermions np(k) is uniform up to the Fermi edge (exclud- 
ing the zero momentum mode at fco = 0), and does not 
depend on the randomly chosen phases 9 m : 



nF,Ds(k m ) = n F ^ DE {k m ) 



if 1 < \m\ < N (24) 
otherwise 



Namely, the SP eigenstates sin(fc m x + 6' m ) are NOs of the 
fermionic system. The Fourier power spectrum of each 
SP state sin(fc m x + 9 m ) [which determine the fermionic 
momentum distribution via Eq. (fTf|) ]. does not depend 
on on the phase 9 m . Each fermionic NO sin(fc m x + 9 m ) 
can be written as a superposition of two plane waves 



sin(fc m x+6' m ) = (e 



ik m x-\-i 



-ik m x—i0 v 



)/2i. Evidently, 



the mean value of the momenta pointing in the positive 



0.8 








-10 k 10 



FIG. 7: Momentum distributions for 10 different states ipDE 
[Eq. (|23|) ] chosen at random, by randomly choosing 10 sets 
{9m | m = 1, . . . , AT} of phases (see text for details). 



(negative) direction is 7r(iV + 1)/L [— n(N + 1)/L, re- 
spectively], that is, it is identical to k pea k- When the 
fermionic states are mapped to the TG states, a wave 
function component which distinctively populates mo- 
mentum modes at ±fc pea fc can appear. This occurs when 
the phases 9 m act coherently, i.e., it is evident that a 
random choice of the phases m destroys the observation 
of the two spikes connected with this component. 

Before closing this section we should say that in all our 
numerical calculations, the phases of the states i/jde were 
chosen at random (with respect to the uniform probabil- 
ity density in [— 7r,7r]). A random choice of the phases 
yields a typical state i/'d e in the sense that one-body ob- 
servables, such as the momentum distribution, of typical 
states approximately coincide. In order to verify this as- 
sumption Figure [7] displays momentum distributions for 
10 eigenstates ipDE (N = 11 particles), obtained by 10 
randomly chosen configurations {8 m | m = 1, . . . , N} of 
the phases. The momentum distribution only slightly 
varies from case to case with one exception that exhibits 
2 (relatively small) dark solitonic spikes. Exceptions from 
the typical behavior will be harder to see for larger values 
of N, because in this case the parameter space spanned 
by N phases 9 m is larger, and it is harder to correlate the 
phases by chance, which could yield characteristic soli- 
tonic spikes in the momentum distribution. Hence, we 
can conclude that our observations regarding the class of 
states ipDE from Eq. (f2"3"]) hold for practically all of these 
states in the sense stated above. 
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FIG. 8: (color online) The RSPDM (a), momentum distri- 
bution (b), and occupancies of the NOs (c) for a DS state 
(iV = 10) in a parity-invariant well-shaped potential, (d) The 
momentum distribution (black dotted line) in comparison to 
the contribution from the 9th and 10th NO: K\$i(k)\ 2 
(red solid line). 



riodic lattices 1131. well-shaped potentials [13(, and so- 
forth. In Ref. 13] it was shown that by parity invari- 
ant filtering of the many-body wave function, one could 
in principle excite the TG gas into a DS state. Let us 
compare the RSPDM and the momentum distribution of 
DS states on the ring, and in a parity invariant poten- 
tial V c {x) = V c a {2 + £ i=1>2 (-) i+1 tanha^Or + (-) l x c )} 
(V c ° = 15, x w — 8, and x c — 25). In such poten- 
tial, DS states are constructed by populating the first 
N odd SP eigenstates on the Fermi side of the map. Fig- 
ure IHa) displays the RSPDM, while Fig. \jgb) displays 
the momentum distribution of such an excited eigenstate 
for N = 10. We clearly observe that the structure of 
the RSPDM and the momentum distribution is similar 
to that of DS states on the ring; the RSPDM has off- 
diagonal mirror-point correlations, while the momentum 
distribution has two sharp spikes. Furthermore, Fig. [5fc) 
displays the occupancies of the NOs, which clearly ex- 
hibit a large and sudden drop after the Nth NO. Fig. &[d) 
shows the contribution from the (N — l)th and the iVth 
NO to the momentum distribution: 53i=g Ai<£»* (fc)<£>^ (fc); 
clearly these NOs are responsible for the peaks in the 
momentum distribution. 



IV. DS STATES IN A PARITY INVARIANT 
WELL-SHAPED POTENTIAL 

The concept of dark " solitons" can be extended to var- 
ious types of parity-invariant potentials (e.g., s ee 113] ). 
DS states are found in harmonic confinement [12j pe- 



The observations presented in Fig. [5] suggest that the 
behavior of the one-body observables of DS states, such 
as the two sharp spikes in the momentum distribution 
and correlation between the mirror points in the RSPDM 
can be found in various types of parity-invariant poten- 
tials. 
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CONNECTION TO INCOHERENT LIGHT 



VI. SUMMARY 



We would like to point out that the behavior of inco- 
herent light in linear [49[ and nonlinear [E(J HJ, [H, [53| 
optical systems has many similarities to the behavior of 
interacting; (partially condensed or non-condensed) Bosc 
gases [m , \53i , l54l . [55j . The dynamics of incoherent light in 
nonlinear systems attracted considerable interest in the 
past decade since the first experiments on incoherent soli- 
tons [50| in noninstantaneous nonlinear media were con- 
ducted. A number of important results were obtained 
(for a review e.g., see Ref. [HT|] ) since then. Among 
the recent results one finds, e.g., the experimental obser- 
vation of incoherent solitons nonlinear photonic lattices 
52 1 , and thermalization of incoherent nonlinear waves 
53j | . We believe that many of the phenomena observed 
with incoherent light in optics [HO, HH, HH, HH can find its 
counterpart in the context of Bose gases. 

In Ref. [l3[ it has been pointed out that there is math- 
ematical relation between the propagation of partially 
spatially incoherent light in linear ID photonic struc- 
tures and quantum dynamics of a TG gas. More specifi- 
cally, the correlation functions describing incoherent non- 
diffracting beams in optics [49] can be mapped [l3[ to 
DS states. However, it should be emphasized that the 
spatial power spectrum of these incoherent beams corre- 
sponds to the momentum distribution of noninteracting 
fermions, i.e., it profoundly differs from the momentum 
distribution of DS states in a TG gas discussed here. 



We have employed a recently obtained formula [40| to 
numerically calculate the RSPDM correlations, natural 
orbitals and their occupancies, and the momentum dis- 
tribution of dark "solitons" in a TG gas. We have found 
that these excited eigenstates of a TG gas have character- 
istic shape of the momentum distribution, which has two 
distinguished sharp spikes; while most of the paper is de- 
voted to the ring geometry, where the spikes are located 
at k P eak — ±7r(A + 1)/L (N is the number of particles 
and L is the length of the ring), we have shown results 
which suggest that such behavior is general for DS states 
in parity invariant potentials. It has been shown that 
the spikes in the momentum distribution are closely con- 
nected to the cross-diagonal oscillatory long-range corre- 
lations between mirror points (x and —x) in the RSPDM. 
This behavior of DS states follows from the fact that they 
are specially tailored; in the ring geometry, it has been 
shown that the two spikes and a special form of spatial 
coherence are lost for most eigenstates that are degener- 
ate to DS states. 
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